from PV_diagram_3 import *

rotation = 190.0 #degrees
Mass = 2.6 #Solar masses 
Mass_err = 0 #Solar masses
v_source = 11.0 #km s^-1
distance_source = 392.0 #pc
inclination = 72.5
ra=83.86514282
dec=-5.15956705

#  ----------------------------------------------------------------
#H2CO
#-------------------------------------------------------------------------------
Object = 'HOPS-370'
Directory = 'HOPS-370/'
#0.25" pixels; 5.25" width
width = 15#pixels	
Line = '$^{13}$CO (2-1)'
File = 'HOPS-370_13CO_image.image.fits'
PV_diagram(ra,dec,File, Directory, rotation, width, Object, Line, inclination, Velocity_Curve = True, mass = Mass, mass_err = Mass_err, v_source = v_source, d_source = distance_source, Thindisk = False, Zoom = True, v_width = 15.0, arcsec_width = 3.0,Overlay_Contour = 'None',imagemin=-99.0,imagemax=-99.0,contour_interval=5.0)



Line = 'H$_2$CO (3$_{0,3}$-2$_{0,2}$)'
width = 15#pixels	
File = 'HOPS-370_H2CO-L_image.image.fits'
PV_diagram(ra,dec,File, Directory, rotation, width, Object, Line, inclination, Velocity_Curve = True, mass = Mass, mass_err = Mass_err, v_source = v_source, d_source = distance_source, Thindisk = False, Zoom = True, v_width = 15.0, arcsec_width = 3.0,Overlay_Contour = 'None',imagemin=-99.0,imagemax=-99.0,contour_interval=5.0)
#0.1" pixels; 2.1" width

File = 'HOPS-370_H2CO-H1_image.image.fits'
Line = 'H$_2$CO (3$_{2,2}$-2$_{2,1}$)'
PV_diagram(ra,dec,File, Directory, rotation, width, Object, Line, inclination, Velocity_Curve = True, mass = Mass, mass_err = Mass_err, v_source = v_source, d_source = distance_source, Thindisk = False, Zoom = True, v_width = 15.0, arcsec_width = 3.0,Overlay_Contour = 'None',imagemin=-99.0,imagemax=-99.0,contour_interval=5.0)
#0.1" pixels; 2.1" width

Line = 'H$_2$CO (3$_{2,1}$-2$_{2,0}$)'
File = 'HOPS-370_H2CO-H2_image.image.fits'
PV_diagram(ra,dec,File, Directory, rotation, width, Object, Line, inclination, Velocity_Curve = True, mass = Mass, mass_err = Mass_err, v_source = v_source, d_source = distance_source, Thindisk = False, Zoom = True, v_width = 15.0, arcsec_width = 3.0,Overlay_Contour = 'None',imagemin=-99.0,imagemax=-99.0,contour_interval=5.0)
#0.1" pixels; 2.1" width

Line = 'SO'
File = 'HOPS-370_SO_image.image.fits'
PV_diagram(ra,dec,File, Directory, rotation, width, Object, Line, inclination, Velocity_Curve = True, mass = Mass, mass_err = Mass_err, v_source = v_source, d_source = distance_source, Thindisk = False, Zoom = True, v_width = 15.0, arcsec_width = 3.0,Overlay_Contour = 'None',imagemin=-99.0,imagemax=-99.0,contour_interval=5.0)
print(v_source)
#0.25" pixels; 5.25" width
Line = 'CH$_3$OH'
File = 'HOPS-370_CH3OH_image.image.fits'
PV_diagram(ra,dec,File, Directory, rotation, width, Object, Line, inclination, Velocity_Curve = True, mass = Mass, mass_err = Mass_err, v_source = v_source, d_source = distance_source, Thindisk = False, Zoom = True, v_width = 15.0, arcsec_width = 3.0,Overlay_Contour = 'None',imagemin=-99.0,imagemax=-99.0,contour_interval=5.0)
print(v_source)

#0.25" pixels; 5.25" width
Line = 'C$^{18}$O (2-1)'
File = 'HOPS-370_C18O_image.image.fits'
PV_diagram(ra,dec,File, Directory, rotation, width, Object, Line, inclination, Velocity_Curve = True, mass = Mass, mass_err = Mass_err, v_source = v_source, d_source = distance_source, Thindisk = False, Zoom = True, v_width = 15.0, arcsec_width = 3.0,Overlay_Contour = 'None',imagemin=-99.0,imagemax=-99.0,contour_interval=5.0)
print(v_source)
#0.25" pixels; 5.25" width

width = 31#pixels	
rotation = 100.0 #degrees
Line = '$^{12}$CO (2-1)'
File = 'HOPS-370_12CO_image.image.fits'
PV_diagram(ra,dec,File, Directory, rotation, width, Object, Line, inclination, Velocity_Curve = False, mass = Mass, mass_err = Mass_err, v_source = v_source, d_source = distance_source, Thindisk = False, Zoom = True, v_width = 150.0, arcsec_width = 20.0,Overlay_Contour = 'None',imagemin=-0.1,imagemax=-99.0,contour_interval=5.0)
print(v_source)
#0.25" pixels; 5.25" width
